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1. Introduction 

This article continues the semiclassical analysis that was started in [I] and |17j . 

In [17, 18J, a general theorem about the construction of solutions of Volterra integral 
equations was proved; [17] also outlined an application to semiclassical approximation for 
the time-dependent Schrodinger equation (converted to an equivalent integral equation). 
The main idea [3j |8] is to use the WKB approximation to the quantum Green function 
(propagator) as the foundation for the solution of the integral equation by iteration. In 
[T?] Theorem 12], which we call the "semiclassical Volterra theorem", it was proved that 
this series solution will converge under two technical hypotheses. However, the issue of 
whether those hypotheses are satisfied in particular concrete problems was not studied 
there. 

In the present paper we apply the semiclassical Volterra theorem to the family of 
potentials z a 8(z) in R 1 and test the validity of the hypotheses. This study was begun 
in [1] , which also introduced the corresponding potentials in R n as a model of boundary 
effects on the vacuum energy in quantum field theory. Here we consider in detail the 
simplest two cases, a = 2 and a = 1. 

The semiclassical construction is based on the classical paths (of the dynamical 
system in question) between two space-time points. Therefore, the first step is the 
classification of these trajectories. For example, if the starting and ending point are 
both in the potential-free region (left of the vertical axis), there is always a direct 
path between them that never enters the potential; its contribution to the semiclassical 
propagator is just the free quantum propagator. For the potentials under study, there 
is usually one other zeroth-order path, which enters the potential region and bounces 
back. We find that this path goes through a caustic (focal point), so that the WKB 
approximation to that path's contribution to the quantum propagator must be modified 
by a Maslov index. Similarly, we analyze the possible paths when one or both of the 
points is inside the potential. 

For each classical path we construct the action and amplitude functions that define 
the contribution of that path to the semiclassical approximation to the propagator. The 
exact propagator is a solution of the time-dependent Schrodinger equation; applying 
the Schrodinger differential operator to any of the semiclassical terms leaves a residual 
expression proportional to the Laplacian of the amplitude function. The proof of the 
semiclassical Volterra theorem shows how to pass from the residual to an estimate on the 
error of a semiclassical solution. Furthermore, when the Laplacian function is bounded 
(no caustics are encountered in the region concerned), a first-order (in K) correction can 
be calculated by concatenating two classical paths and integrating a certain integrand 
over the location of the point where they are joined. If the potential is smooth, this 
process can be continued to arbitrarily high order, although the formulas rapidly become 
very cumbersome, and the theorem shows that this series in h is convergent to the exact 
solution. 

In the model under study this idealized strategy is eventually disrupted in two 



WKB Approximation to the Power Wall 



3 



ways. First, our model potential has a mild singularity at z = 0, so the semiclassical 
expansion fails at a sufficiently high order. Second, and more fundamentally, caustics 
do occur, even where the potential is smooth, and they cause blowups in the residuals, 
error estimates, and higher-order corrections. The ultimate remedy for this disease is to 
convert the calculations to momentum space [E2, HU [26], a project that goes far beyond 
the scope of the present paper. 

Section [2] reviews general semiclassical (or WKB) theory, and section [3] sets up a 
class of models, the power walls, and begins the analysis of their classical paths. Detailed 
semiclassical treatments of quadratic and linear power walls are presented in section H] 
and section [5J respectively, and section |6] presents some conclusions. 

2. Hamilton— Jacobi Theory, WKB Approximation and Volterra Integral 
Equation 

In this section we review the general theory of construction of the semiclassical 
propagator A sc i(x, t; y, s), following such treatises as [6] and [T5] . 

Consider a quantum particle subject to a (sufficiently smooth) potential V(x, t), 
x G R™. A natural ansatz for the wave function is 

V>(x,*) = A(x,t)e^ (x ' 4) (1) 

where A(x, t) and S(x, t) are called the amplitude and the action of t), respectively. 
Substituting (CQ) into the time- dependent Schrodinger equation 

one obtains the partial differential equation 



= A 



dS 1 



dt 2m 



dA 1 



+ — (vsy + v -m — + — (va- vs) + — aas aa, (3) 



. dt m 2m 



■i 



2m 



where A is the Laplacian operator. Separate the real and imaginary parts of ((2]) to get 
dS 1 ,„ m2 T7 h 2 AA 

m + 2^ (VS > ^ 2m~A~ < 4 > 

and 

BA 1 
m— + (VA ■ VS) + -AAS = 0. (5) 
dt 2 

The classical limit is obtained by taking the limit h — > 0, whereupon (j4j) becomes 

| + J_ ( v S )S + V(x, t )=0. (6) 

This is the Hamilton- Jacobi equation. The phase S"(x, t) is interpreted as the classical 
action. 

Equation ((6]) has the form || + H(x,VS(x,t),t) = 0, where H is the classical 
Hamiltonian function, H (x, p, t) = ^\p\ 2 + V(x., t). A classical solution is a local curve 
x(t) satisfying the equations 

dx(t) dH p 1 



dt dp m m 



V5(x(t),t) (7) 
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and 

Equation (CO) enables one to construct the action S(x, y, t) from a knowledge of a family 
of classical solutions x(t). The total time derivative of the action must be 

^=^+x-VS = -tf + x- p = L(x(t), x(t)). (9) 

This equation implies that we can get solutions of the Hamilton-Jacobi equation by 
integrating the Lagrangian L along the trajectories: 

S(x,y,t) = / L(x(u), x(u)) du + S , (10) 



Jo 

where So is initial data, and S(x,y,t) then solves the Hamilton-Jacobi equation. 

Conversely, if we have a local solution of the Hamilton-Jacobi equation, then 
VS = p = mx along its classical trajectories. Then to solve the Schrodinger equation 
through order h we still need to solve the transport equation ([5]). The latter can be 
rewritten as 

-±AAS=(^ + ±VS.v)A=(i+±.v)A= d 4, (11) 
2m \dt m J \dt J dt y ' 



or 



1 1 dA d , . , 
AS = -— = — In A. 12 

2m A dt dt K ' 

We can solve for In A by integrating along the classical trajectories, with the result 



1 



A(x., t) oc. exp / AS(x(m), u) du (13) 

L 2m J J 

(where the y dependence has been suppressed). However, the amplitude function can 
be expressed in an alternative way: 



4(x, t) oc Vdet C, (14) 

where C = V x V y S ; detC is known as the Van Vleck determinant. The fact that this 
determinant is a solution of the transport equation is well known but nontrivial [6]. 

In what follows we restrict attention to time-independent potentials V(x) and 
concentrate on the quantum propagator, the Green function that produces a solution of 
the time-dependent Schrodinger equation from arbitrary initial data. K can be written 
in Dirac notation as 

K(x,t;y,s) = (x.\U(t,s)\y), (15) 

where U(t, s) is the unitary time-evolution operator for the system taking states at 
time s to states at time t . Intuitively, the propagator i^(x, t; y, s) describes the motion 
of a quantum-mechanical particle travelling from the space-time point (y, s) to the 
point (x, t) and can be interpreted as passing through each possible intermediate point 
(r, r) with a certain probability amplitude. The basic concept that underlies the theory 
of higher-order semiclassical approximations is that in a local space-time region the 
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particle evolves under the semiclassical propagator between encounters with the effective 
potential V sc i = ^r- (This idea was developed by Balian and Bloch [3] in the context of 
the time-independent Schrodinger equation and its Green function, and mentioned by 
Dowker [8j sec. 3.2] for the time-dependent case.) 

For time-independent V, K is a function of the time difference t — s, so one usually 
sets s = without loss of generality, and the fourth argument of K is suppressed in the 
notation. (Similarly, the action and amplitude functions in full generality are functions 
of s as well as of (x, t and y.) The Green function if (x, t; y, 0) satisfies the homogeneous 
Schrodinger equation in the variables (x, t), except at the source point (y, 0). Therefore, 
the machinery introduced above applies to it. By a standard argument (e.g., [11]). it 
can also be defined by the nonhomogeneous equation 

d\ 

H x - ih— if (x, t; y, s) = -zM(x - y)S(t - s), (16) 



9t. 

where H x is the Hamiltonian of the quantum system, appearing here as a function of 
the x variable, and S is the Dirac delta function. Under suitable technical conditions 
the solution of the nonhomogeneous equation 

H + ih^ji}(x,t) = <t>(x,t) (17) 

with initial data ^(x, 0) = is 

V(x,f)=/ / K{^t-y,T)4>{y,T)dydT=[K<j)}{^t). (18) 

JO JR n 



With respect to t, the integral operator in (118]) is of the Volterra type; that is, the upper 
limit is the solution's variable t, rather than +oo or some large fixed T. 

The free propagator iff (x, t; y, s) = iff (x, y, t — s) in the space-time M n x M + with 
V(x, t) = is well known to be (for t > s) 

/ \ n/2 

K t (x,P,y,s)= U-^—,) e H*-yP/»(i-.>. (19) 

The exponent in if/(x, y, t) is j times the action ^(x, y, t) for a free particle. The Van 
Vleck determinant for this case is 

detf-^M = (20, 



dxidyj/ V t 

Thus the free quantum propagator fits into the WKB framework as 



/ 1 x n/2 / , f)2 C s f A >. 

K l{ x,y, t) = Mt )e^=(— ) ^ det (- 5 -J-)exp{ s5o (x, y , t )}. (21) 

The particular normalization factor in (fT9|) is the one that gives the correct initial value 
to if(x, t;y, s) on the surface t = s. Alternatively, if one thinks of K(x,t;y,s) as a 
solution of the nonhomogeneous Schrodinger equation (Tl6]) in all of space-time, it gives 
the correct delta-function singularity at (x,t) = (y,s). 
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Similarly, we will be able to write a semiclassical propagator K sc i(x., t; y, s) in the 

form 

K sc i(x, t- y, s) = (2nifr)- n/2 VdrtCe lS/h (22) 

where C is the nx n matrix with elements Cjj = — d ® 2 g y . The factor V det Ce lS ^ h arises 
as the solution of the transport equation (J5J), and the arguments for the normalization 
factor are the same as in the free case. 

A well-known technique [T0| [T3| [22] for the construction of Green functions for the 
Laplace and Helmholtz equations, and also the heat equation, in bounded domains in M. n 
(billiards) is by reduction to integral equations on the boundary. An important feature 
of the heat equation is that the solution of the boundary integral equation by iteration is 
convergent because of its Volterra structure. Therefore one has, in principle, an explicit 
construction of the solution. The Schrodinger equation has the same Volterra structure, 
so one expects again to have a convergent series solution. A general theorem to this 
effect was proved in [TTJ, Ch. 6] with corrections in [18]. The application of this general 
Volterra theorem in a particular context reduces to showing that the operator family 
arising in the construction of the individual terms in the series is uniformly bounded on 
a suitable Banach space. In any particular case this may be a nontrivial task and may 
require additional technical assumptions. 

In the semiclassical Schrodinger problem fl2]), the key idea that we implement is 
to use the WKB approximation to the quantum kernel analogously to the free kernel 
approximation in billiard problems. In the billiard problem, scattering happens only 
at the boundary; in our case, the particle is scattered throughout the bulk region by a 
source that is the residual error in the WKB approximation to the exact kernel. This 
construction is developed in [TTJ Ch. 8]. The WKB kernel is 

K scl (x, t; y, 0) = (2mh)- n / 2 Ae' s/h , (23) 

where, as explained above, 

S(x, y, t) = L(q(r), q(r)) dr, L(q(r), q(r)) = yq 2 - V(q), (24) 
is the classical action, and the amplitude A is 



d 2 S 

• 4(x - y - () = '/ det ^&-- (25) 

If there is more than one classical trajectory q(r) starting at y at time and 
arriving at x at time t, the semiclassical approximation is a sum of such terms, possibly 
modified by Maslov phase factors to keep track of places where the radicand in (1251) 
becomes negative. 

We define a kernel Q by 

Q(x, t; y, r) = (2mhy n / 2 H 2 [A x A(^ t; y, T )}e lS ^^' h . (26) 
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Both Q and K Bcl define Volterra operators by formulas precisely analogous to (fTHl) ; we 
denote these operators (including both time and space integrations) by the same letters 
as the kernels. The operators K, Q and K sc \ are related by 

K- X K^ = (-iM t - h 2 V 2 + V(x))K 8cl = 1 - Q ; (27) 

that is, Q = 0(h 2 ) is the amount by which K sc \ fails to solve the PDE for which it was 
devised. Thus, formally 



(28) 



k = k sc1 (i-q)- 1 = k sc1 J2Q 3 - 

3=0 

The Volterra property ought to enable one to prove that this series converges. 

To examine this question in detail, we write out the kernel representation of (J25 
We introduce a notation for the spatial parts of the integral operators: 



[r(t,r)0](x)^[^ scl (t,r)^](x) 



K Bc i(-x,t;y,r)(f)(y,T) dy, 



[Q(t,r)0](x) 



<2(x,t;y,r)0(y,r) dy. 



Note that Q from (|26|) can be written 
Q(x,t;y,r) = h 



K sc i(x, t; y, r) 



so 



[Q(t,r)0](x) = /l 2 



A x A(x,t;y,r) 



r A(x,t;y,r 
Let ip = K<p and apply (T2"8l . The zeroth-order term is 

JW(x,t)= / [r(t,r)0](x)dr. 
Jo 



-^sci(x,t;y,r)0(y,r) dy. 



(29) 
(30) 

(31) 
(32) 

(33) 



Note that if T(t,r) is a bounded operator on L 2 (R"), with bound C independent of t 
and r, then 

||-Kbc10||l3(R»)(*) < C / H0l|L 2 (R™)( r ) ^ r < iC||0lU°°.2(7 t ,R") , (34) 

Jo 

where the final norm is the supremum of ||</>||L 2 (R™)(' r ) over t E It = [0, t}. (For simplicity 
of notation we consider only t positive.) The first nontrivial term in the series is 



lf sc iQ0(x, t) 



dnT(t,Ti) 



dr[Q(ri,r) 



UO 



h 2 / dri / dxiKgdfat^iiT!) x 
A Xl A(x 1 ,r 1 ;y,r 



dr dy 
o V 4(xi,t;y,r) 



^ S ci(xi,t;y,r)0(y,r). 



Now suppose that 

A x A(x,t;y,r) 
^(x,t;y,r) 



< .D for all x, t, y, r. 



(35) 



(36) 
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Then 



\K scl Q<p{x,t)\ < h 2 D 



dn cbni^x^x^n) x 

JR" 

dr / ciyi ; r sc i(xi,t;y,r)0(y,r) 

h 2 D f d n r i c /r[r(t,r 1 )r(t 1) r)0](x) 

Jo Jo 



and so 

||# sc iQ0|U W (*) < ^ 2 ^C 2 / dn / 1 dr ||0|| L2(R ™(r) < ^DC 2 ^ . (37) 

Jo Jo ^ 

Similarly, or by induction, for each j 

ll^sciQ J ~VlU w 0O < tFiy-W | u\\ L ^ {Ium . (38) 

The sum of all such terms is majorized by an exponential series, so it converges in the 
topology of L 00,2 (Jy, M n ) for any oo > T > t. We have thus proved the following. 

Theorem 1 Suppose that the following two hypotheses hold: 
AA 
~A~ 

(ii) K sc \(t,r) (defined in (dSjjj is a uniformly bounded operator from L 2 (R") to itself. 



(i) 



< oo (cf. (Esp;. 

L°°(/r 2 ;K 2n ) 



Then for all <fi G L°° ,2 (It, ffi") the inequalities (E^j hold. It follows that the series 



oo 



K<f ) = Y,K sc iQ 3 <p (39) 

converges (for t < T). In other words, the nonhomogeneous time- dependent 
Schrddinger equation (T7\ ) can be solved by iteration starting from the semiclassical 
kernel approximation, K sc \ f j2B(l . (Efy), to the full propagator, K (see 

This theorem was proved in [TTl Chapter 8] . The idea of proof is the same as that for 
a general theorem on the solution of Volterra integral equations by iteration, formulated 
in [TTl Chapter 6] with a flawed proof and successfully proved in [18]. Because of the 
complicated structure of the expressions K sc \Qi , it has been more convenient to repeat 
the argument from the beginning rather than to force the problem into the mold of the 
general Volterra theorem. 

Because Q(x, t; y, r) for fixed y is not an L? function of x, we cannot apply the 
argument to get literal convergence of the series for -^(x, t; y, r) (without the "smearing 
function" 0). Moreover, the convergence is in L 2 , not pointwise. Of course, it is possible 
that a stronger theorem holds. 

The condition (ii) is expected to hold wherever K sc \ is a decent approximation to K 
(i.e., away from caustics). The point is that for fixed times the operator defined by K is 
unitary (cf. ( jT5l) ). so that the one defined by K sc \ should be approximately unitary and 
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hence bounded. This reasoning would not apply if L? were replaced by L°°; indeed, the 
propagators are not bounded functions as t — > 0. 

The determinant in ( 125]) is singular at caustics, where the mapping from initial 
velocity data (at y) to x ceases to be a diffeomorphism. One can expect both conditions 
in Theorem [1] to become problematic if the orbit goes through a caustic. A way to 
go beyond caustics (if necessary) is provided by the Maslov theory [15], as already 
implemented in a similar problem in [24] . 

The construction in Theorem [T] implements the Feynman path integral idea in a 
way different from the usual time-slicing approach. (A similar observation was made by 
Putrov [21] in a different context.) Each term in ( 1391) is an integral over classical paths 
with j scatterings off an effective potential A A/ A. 

3. The Power Wall Potential 

3.1. The Model 

To test the validity of the two hypotheses in Theorem [1] in the context of a concrete 
problem, we consider a family of potentials in one dimension, namely 



where 1 < a G K.. The study of this "power wall" model, in a wave equation, was 
initiated in jl] and continued in [19] in the context of quantum vacuum energy, and we 
hope that our study of the associated Schrodinger equation will yield new information 
about the spectral density (and hence the vacuum energy) of the operator —V 2 + V . 
(In the vacuum-energy papers there were two transverse dimensions, but here we ignore 
them because their contribution to the quantum kernel in dimension 3 is a trivial factor.) 
In [1] the coupling constant A was written for any a in terms of a dimensionless constant 
and a fundamental length, but here we are concerned with particular values of a and 
will choose the physically most natural notation in each case. 

The most calculationally tractable values of a are 2 and 1. They are investigated 
in detail in the next two sections. 

3.2. Classification of the Classical Paths 

A particle moving in the potential PO]) is acted on by a force that never points to 
the right. Therefore, the possible trajectories have only a small number of possible 
"topologies". Consider first the initial-value problem, where q(0) = y and q(0) = v = 2p 
are prescribed and one solves for q(r). If y < 0, the particle initially moves freely; if 
also p < 0, it will move freely forever, but if p > it will eventually enter the region 
with the potential. In the latter case it will accelerate to the left and eventually exit 
from the potential and move freely again. If y > 0, the particle immediately accelerates 
to the left and eventually reaches the free region. 




if x > 0. 



if x < 



(40) 
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r 












y 



q 



Figure 1. Classification of paths with y > 0. Solid: Type D. Heavy: Type B. 




y 

Figure 2. Classification of paths with y < 0. Dashed: Type A. Solid: Type E. Heavy: 
Type C. 

In semiclassical analysis one needs to solve the two-point boundary-value problem 
where q(0) = y and q(t) = x are prescribed. It follows from the foregoing remarks that 
five types of path are possible, which we letter in order of increasing complexity and 
indicate in Figures HH2J 

• For y > 0: 

If x > 0: Type B 
If x < 0: Type D 

• For y < 0: 

If x > 0: Type C 

If x < 0: Type A or Type E 

As the last condition shows, it is possible for two points in space-time to be joined 
by more than one classical path. Conversely, for certain values of y, x, and t it is possible 
that a path of a certain expected type will not exist. The precise constraints on the 
parameters will depend on a. (For example, when a = 2 a complete excursion within 
the harmonic oscillator potential must take an elapsed time of precisely half a period, 
thus in type E there is a lower bound on t that does not apply when a = 1.) 

For more refined semiclassical approximations it will be necessary to solve the two- 
point boundary-value problem with p, x, and t given. For each such data list, the sign 
of y may not be immediately obvious, so one must explore several of the five types for 
existence of paths. This problem is left for later work. 

All the paths discussed so far are "zeroth-order" paths, needed to construct the basic 
WKB propagator i\~ sc i . The first-order approximation ( |35i) will introduce concatenations 
of paths of types C and D, with a possible change of velocity at the joint, whose 
contributions must be integrated over the location of the joint. Furthermore, because 
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our model potentials are not smooth at x = 0, there is another sort of path that must 
be considered in a complete semiclassical treatment. These are the paths that reflect 
off the singularity at the origin. The contributions of such terms decrease rapidly with 
a (at least for integer a), but for a = 1 (considered in section [5]) they are at least 
comparable to ( 1351) . We hope to return to them in later work. 



4. The Harmonic Oscillator and the Quadratic Wall 

We now consider in detail the model with a = 2, 

( if x < 0, 

V V) = 1 2 2 . f > n ( 41 ) 

-uj x it x > 0. 
v 4 

To simplify the formulas we take the mass m to be |. As previously noted, each classical 
path contributes to the leading-order semiclassical approximation to the propagator. 
We calculate the action S, the amplitude A, and the Laplacian AA. The Laplacian 
is of interest because it is a crucial factor in the source term for the next-order 
approximation (see f )26|) and f )35|) ). In other words, A A/ A is the residual in the leading- 
order approximation (the right-hand side of (HI)); a singularity in it, in particular, signals 
a breakdown in the approximation. 

Because we are considering a system with time-translation invariance, as explained 
in section [2] we usually use the notation K(x, y,t — r) for the quantum propagator 
K(x,t;y,r). But we must retain the initial time variable r to facilitate concatenation 
of paths. 



4-1- Type A: Free Particle 

Whenever x and y are both negative (i.e., the initial and final points are on the left side 
of the origin, Figure [3]), there always exists a direct path between them that stays in 
the force-free region. The particle's position and velocity then are 

ff(r)=» + r(^), v(r) = q(r) = . (42) 

The classical action and amplitude of the free direct path are thus 
n f*l. 2 , {x-yf a2 d 2 S 1 

S= -q 2 dr = v y , A 2 = -— — = — . (43) 



4 At dydx It 

(We note that the Lagrangian for a path of this type is equal to the (constant) kinetic 
energy of the particle.) The WKB construction ( 123]) therefore yields the propagator 
(fl9l) (with n — 1 and m — |). The Laplacian A A is 0, because in this case the WKB 
propagator is exact: 

-| ^+00 /TON 1/2 

K f {x,y;t) = dke ^-y) e -^t/(2 m) = e — _ (44) 




4-2. Type B: Harmonic Oscillator 

Now let x and y both be in the potential region (i.e., to the right of the origin, Figure 
H]). Eventually, the boundary cases where one of the coordinates is equal to will be 
included. 

The general solution for a particle that remains inside the potential is 

q{r) = acos(wr) + frsin(wr). (45) 

Consider the initial condition 

V = 9(0) = a (46) 

and the final condition 

x = q(t) = ycos(ut) + 6sin(a;t), (47) 

which implies 

x — y cos (cut) 



sm(ut) 

For future use we restate the solution for a general starting time, s: 



( 4 * 



i \ , , x — ycos(co(t — s)) . , , lAnS 

q(r) =ycos(u(T - s)) + y - V V ;; sm (w r - s . 49 

sin(w(t — s)) 



(We always assume s < t.) 
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Returning to s = 0, we find the Lagrangian 



,2 



L{x,y,T) i 



x — 2xy cos(ojt) + y cos(2ojt) 



cos(2wr) 



, x — y cos(ut) . 
- 2y U — J sm 2wr 



sin 2 (tut) 

(50) 



and hence the action 

S{x,y,t) = I L(x,y,T)dr 



= — — r [x 2 cos( cot) + y 2 cos(wt) — 2xy] . (51) 

4sin(wt) 

The amplitude is given by 

A 2 = -^-= " x . (52) 
oxoy 2 sin (cut) 

As in the free case, A is independent of x and hence AA = 0. So the WKB propagator 
for this type of path is exact and equals 



K no (x,y,t) =,/ - , — exp{- — [x 2 cos(ut) + y 2 cos(ut) - 2xy] }. (53) 
V 47TZ sin(o;t) I 4j sin J 



47rz sin(a;t) I Ai sin(a;t) 

This is the well known [HI [151 El] quantum propagator for the one-dimensional harmonic 
oscillator for < t < tz/lo. 

We pause now to note the consequences of the fact that the period of a harmonic 
oscillator is independent of the amplitude. A trajectory of the sort shown in Figure 
Upright), starting and ending at the origin, exists only when t — s = tc/u, and in that 
case it is not unique (the amplitude is arbitrary). If either x or y is strictly positive, the 
(unique) path (j49j) is a segment of one of those paths, and thus necessarily t — s < ir/u. 
The formulas ()48l)-(l5TT) remain meaningful and correct when one, but not both, of the 
endpoints is 0. For example, when y = the action just simplifies to 
s( ) = 

4sm(u;t) 

In the special case when x = y = and t = tt/uj, one must return to the general 
solution (145 j) and impose the boundary conditions 

= y = q(0), = x = q(t) (55) 

and a supplementary condition, 

q(0) = v. (56) 

One finds the solution q(r) = bsm(ur) along with the relation b = v/uj , the consistency 
condition t = tt/u, and the final velocity q(t) = —v, which is obvious from conservation 

of energy. In terms of a general starting time, 

v 

q(r) = — sin(a;(r — s)). (57) 
In this special case the Lagrangian is 

L = ^-[foVcosVr) - b 2 oj 2 sm 2 (ojr)}, (58) 
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and so the action is 

If 1 1 

S = -v 2 / cos(2c<;t) dr = —v 2 sin(2u>£); 
4 J 8u 

but t = n/u), so 

S = (59) 

for these special paths. 

The nonuniqueness of the special paths and the accompanying singularity in A at 
t = ii/co can be regarded as a caustic. However, it is a very unusual kind of caustic, 
inasmuch as it does not constitute a breakdown of the WKB approximation. The 
function K^q is an exact solution of the Schrodinger equation when extended through 
the singular point with the correct Maslov phase J2U [15] . 

4.3. Type E 

The next simplest case to analyze is when x and y are both negative but the path goes 
through the potential region (Figure [5]). Such a trajectory is effectively a concatenation 
of three types: one part of the path is the special case of type B and the other two 
parts are instances of type A. Let t\ be the time when the particle first crosses into 
the potential region and £2 be the time where it exits the potential region. Because the 
potential is quadratic, we know immediately that t > - and £2 = £1 + -• Also, 

v (tj) = - V- > and v(t 2 ) = < 0. (60) 

t\ t — ti 

By conservation of energy, v(t\) = — v (£2) and hence 



t 



t 



(h + - 



Therefore, 



V <t-*-\, (61) 



and hence 



x + y V uj 
y(t — — ) 7r 

t 2 = zl <al + _ . 62 

x + y uj 

The trajectory can now be determined from (|5"T1) and the initial datum 



to be 



= - t = - 1 — — ( 63 ) 
9(T) = 7 ^ ± 4sin(u;(£-£ 1 )). (64) 

[Tl — Qjt) 

Now, we compute the action by integrating the Lagrangian along the entire path: 

S E {x, y, t) = S A i(x, y, £) + S B {x, y, t) + S A 2(x, y, t) 

rti pt2 rtz 

= / L(q,q)dr+ / L(q,q,)dr+ / L(q } q) } dr. (65) 
Jo Jti Jti 
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t 


x 


V 





Figure 5. Type E 



But Sb = according to fl59|) . and the other two terms are cases of (}4"5]) . Thus the total 
action is 

S E (x,y,t) = ^ + -7^—, = ^ + V ll ■ (66) 
The amplitude function is then 

A - dxdy ~ 2(t - I) ' {b7) 

and once again AA vanishes, Therefore, the solution for type E is also exact (but this 
string of good luck is about to end). 

The negative sign in f )67|) should not be a surprise. By continuity from the (purely 
harmonic) case (y = 0), when y < but y is small one would expect a caustic to occur 
somewhere near x = —y, t = tt/oj. Therefore, when the trajectory reemerges from the 
potential, this term of the kernel carries a Maslov phase factor of — i [151 ttH 121] • The 
existence of the caustic will be verified in the next subsection. 

44. Type C 

Here we have a concatenation of type A with type B (see Figure El left). Again we let 
ti denote the time when the particle passes from the free region to the potential region. 
For the first segment of the path, the velocity is 

q(t 1 )=q(0)) = -y/t 1 . (68) 

For the segment of the path inside the potential, we use the solution ( I49p found for 
type B, with y — and s = t 1 : 

q(r) = - * sm(u(r - h)). (69) 

sin(w(r — ti) 

Alternatively, since ( 157]) with v given in f )68|) is independent of x, it applies to our path: 
q(r) = -^-sm(uj(T-t 1 )). (70) 

Combining these two equations (or simply setting x = q(t) in the second one) yields the 
relation 

uxt\ + ysm(u)(t — ti)) — (71) 
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Figure 6. Type C (left figure) and Type D (right figure) 



to determine t± . 

The transcendental equation ( T7TT) for t\ cannot be solved in closed form, but by 
rescaling the physical variables to dimensionless quantities, we can get simple master 
plots of the linear and sinusoidal terms, as shown in Figure 7. Intersections of the lines 
with the sinusoid represent solutions to f lTTj) and yield valid classical paths. Let 

6 = uti , p = —y/x, T = ut. (72) 

Then (I7T1) becomes 

(73) 



= psin(T- 0). 
Setting Q = T — Q = u(t — ti) improves the equation further to 
T — Q = psinfi. 



(74) 



Note that t\ < t < ti + - and hence 0<T<0 + 7r, or < f2 < 7r. Thus our task is to 
find zeros, in that interval, of 



/(fi) = psin(Q) + tt -T. 



(75) 



It is clear from the figure that the number of such solutions can be 2, 1, or 0. Observe 
that if 



then 



= f'(Q) = pcos(fi) + 1, 



Q = cos 1 ( 

v P 



and that this situation can occur only if p > 1. The condition for the sine curve and 
the diagonal to be tangent is /(fi) = = f'(Q), whence 

1 

/)// ^ \ P y 



= /(cc. S -(- I)) = psin(cos-(- I)) + cos-(- i) - T 



T. 



Therefore, we define 



V - 1 + cos ^ (- -) 0>1) 



(76) 



and conclude 
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7T 



T 




jp^- 


1 


i 

p 







(b) 



Figure 7. (a) The four possible relations between a diagonal line and the principal arc 
of the sine curve, (b) Resulting division of the parameter plane, labeled by intersection 
numbers. 



T > T* =>- no solutions; 

T = T* =>■ 1 solution; 

T < T* =>- 2 solutions, if T > 7r; 

T < 7r =>- 1 solution. 



The action for a type C trajectory is the sum of the actions for its two segments, 

(77) 



y ujx 2 cos(ut) 

Sc(x, y, t) = — h ... 

4ti 4 sin(wt) 



+ 



sin(2co>(t — ti)), 



(7f 



4tj ' 8wt^ 

where the first version comes from (|54|) and the second follows by (ITTj) . The advantage 
of the second form is that it is independent of x and hence usable at every point on the 
curve. 

We turn to the calculation of the amplitude and its Laplacian. By implicit 
differentiation of ( ITT]) in the form y sin Q + xuti = we find 
dt\ ti dti x ti 

dx 



(79) 



x — y cos Q dy y x — y cos Q 

Computing derivatives of ( 1751) directly leads to complications, but from ( JTTj) we get the 
simple formula 

dS y , , 

dy 2t x 1 ; 

(This result is recognized as the negative of the initial momentum of the particle, as the 
Hamilton-Jacobi theory dictates. It is determined by the initial (free) segment of the 
trajectory regardless of what is considered to be the final endpoint, x; this explains why 
only the first term of (175|) can contribute.) We now easily get 

A 2 = d 2 S _ y dh _ y 1 

dxdy 2t\ dx 2t x x — y cos Q 

There will be a caustic if the denominator of ( 1811) changes sign. (Since that divisor 
arises from dti/dx, its vanishing says that t\ (hence y) can vary without changing x (at 
least to first order). The caustic thus represents a kind of nonuniqueness or degeneracy 
of the family of paths.) To study this issue we consider a fixed trajectory with a moving 
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endpoint (that is, fix y and t\ and let x and t vary). When t « t% , f'(£l) is large and 
positive (p — > oo, cos(f2) — > 1). Near the exit point, t ~ is > /'(^) is large and negative 
(p — >■ oo, cos(fi) — >■ —1). Therefore, every trajectory does pass through a solution of 

= f'(Q) = pcos(fi) + 1, (82) 

which is a singularity of (|8T|) . somewhere on its retreat from the potential. This verifies 
the last claim in section 14.31 

In general it is convenient to find the Laplacian of A from A 2 by defining B = A 2 
and noting that 

VA = ~B~ 1/2 VB, (83) 
2 



so that 



A! -B- 1 AB--B- 2 (VB) 2 . (84) 



A 2 4 
In the present case, flHT]) and f )84l yield 



and 



where 



2(x-|/cos(a;(t-t 1 ))) 3 1 J 

— = ^[4^cosfi - 6x 2 + 2y 2 + 3x 2 u; 2 t 2 ], (86) 



K = x -ycosiuit-tx)), fi = w(t-ti). (87) 

We observe that Y vanishes at the caustic but the numerator of f )86|) does not. Therefore, 
unlike the case ( |53l) . there is a genuine breakdown of the semiclassical propagator 
approximation in the vicinity of the caustic. On the far side of the caustic (larger t) 
the approximation will again be good if A (now imaginary) is assigned the phase —i, in 
keeping with the general Maslov theory [T5"| [T4"] . 

4-5. Type D: Type C Reversed 

A type D trajectory has the form q(r) = a sin(a;r) +b cos(ut) inside the potential. With 
q(0) = y and q{ti) = 0, we arrive at 

q{r) = — y cot(uti) sin(o;r) + ycos(wr). (88) 



The total action of the classical particle, in analogy to ( 1771) . is the sum of that for a B 
segment and an A segment: 

1 x 2 
S D (x, y, t) = -y 2 u cot(uti) + _ — . (89) 

Applying the velocity condition to get an implicit equation for t\ yields 

x 

Q{ti) = 7 \ = —yu cos(uti) cot(wti) — yusm(ujti), (90) 

(t - ti) 
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which leads quickly to 

yuj{t-tx) +xsm(uti) = 0. (91) 

It follows that 

dh = t - ti dh = _y t - ti 

dy y — x cos(ujti) ' dx xy — xcos(uti) 
To find the amplitude function of the trajectory we must take the partial derivative 
of t\ with respect to x. Because |^ is the final momentum, we can obtain the amplitude 
formula from just the second term of So if we take the x derivative first: 

x dti 



A2 d 2 S d 
A 2 = 



x 



dydx dyV2{t-t 1 )\ 2(t-t 1 ) 2 dy 

x 1 , . 

= - r ; 7 ■ 93 

2(t — ti) y — x cos(a;ti) 

From ( |93l) we obtain 

VB = y^tlM (94 ) 

2{y-xcos{ut 1 f ' 1 ' 

and hence by fl84l) 

^ = - y2u 2^ lf [Axyc °* {uti) - 6x2 + 2y2 + yW{t - ti)2] - (95) 

Here Q = u(t — ti) and 

X = y — x cos(utfi). (96) 

The question of when paths of type D exist is very much like that for type C, with 
the roles of x and y interchanged. The lack of symmetry between fl86l) and fl95|) arises 
because the calculation takes place at the endpoint inside the potential in one case but 
outside in the other. 

In future work we hope to tackle the next (single- reflection) term, ( |35|) . in the 
series (1391) by concatenating trajectories of types C and D. The final point of the first 
trajectory is the initial point of the second, but the momenta need not match up, since 
the particle is scattering off the effective potential ( 186|) . For given (x,y,t) outside the 
potential (HTj) . one must integrate over all {q,r) inside the potential for which such a 
path exists. From the taxonomy of paths explained in connection with Figure [7] it is 
clear that as many as four trajectories can exist. Thus the kernel of K SC \Q turns out to 
be a sum of four terms, each with a domain of integration that is a nontrivial subset of 
the region (0, oo) x (0, t). 

More precisely, in the C trajectory let us rename x as q and t as r and add on a 
new D trajectory from (q, r) to (x,t) (with x < and t > r). The role of y is now 
played by x, and that of r (which was formerly t) is now played by t — r. For the 
new trajectory we have a new time parameter T and a new scaled position parameter 
p = —-. There are no solutions if T ^> max{7r,p} or T > max{7r,p}. If T ^> ^, then 
p = — | and p——~. If t 3> ^, these conditions are both satisfied only for very small q 
(large p and p). For small p there are also solutions with large q. The boundary curves 
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Table 1. 


Important boundaries 






Old path 


New path 


P 


y = -Q 


x = — q 


T = TT 


UJ 


t-T= 2L 


T = T* 


r=^V» a -9 a + «»- 1 (f) 


i- T =J-^2_ ?2+cos -l ( i ) 



separating the integration domains are indicated in Table [U Further analysis is deferred 
to future papers. 



5. The Linear Potential Wall 



In this section we consider the wall potential with a = 1: 

V( { ° if X < °' (97) 
\kx if x > 0. 

This problem is in some ways harder than that of the quadratic potential, because the 
period of the motion in type B is no longer independent of the amplitude. Also, we will 
find that the WKB construction is no longer exact for type E. On the other hand, the 
equations that implicitly determine times when the trajectories cross the vertical (time) 
axis are not trigonometric but cubic and hence can be explicitly solved. 

As usual, we set h~ = 1 and m — h, and the variable x will be replaced by q(r) when 
we are calculating an entire trajectory starting at y and ending at x. The equation of 
motion, inside the potential region (positive g-axis), is 

q(r) = -- = -2k, (98) 
m 

with the general solution 

q(r) = -kr 2 + ar + b. (99) 



5.1. Types A and B 

Type A for the linear potential is, of course, exactly the same as before, so the formulas 
in section 14.11 still apply. 

For type B we must use (1991 with prescribed endpoints, 

g (0) = y > 0, q(t) =x>0. (100) 

The results are 

(x — y) + kt 2 

b = y, a= K - U -j . (101) 

but it is convenient to refrain from subtituting the cumbersome expression for a until 
after the following calculations are finished. We have 

q(r) = -2kT + a (102) 
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and hence 

q 2 = 4k 2 T 2 - Akra + a 2 , V(q) = -k 2 T 2 + kar + ky. 
Thus the Langrangian is 



L(x,y,r) = -q 2 - kq 



1 



2k 2 r 2 - 2kra - ky + -a 2 , 



(103) 



and the action for this path is 



S(x,y,t) = / L(q,q)dr 
Jo 

*2k t t 2 ■ 

= kt a H kyt, 

3 A y ' 



a 2 t 



which becomes after a final simplification 

(a; — y) 2 (x + y)/ct /c 2 t 3 
~ it 2 12~ ' 

Only the first term gives a nonvanishing contribution to 

d 2 S 1 



A 1 



dy dx It 

so the amplitude is the same as for a free particle, and A A = 0. 
It follows that the WKB propagator for this type is 



(104) 



(105) 



K L (x,y,t) 



Anit 



exp 



rz(x - y) 2 - 




" i(x + y)kt~ 




[ ik 2 t 3 ^ 


V At . 


exp 


i 2 J 


exp 


i 12 J 



(106) 



and that it is exact. Indeed, fl 1 6 [) is a well known expression for the quantum propagator 
of a particle in a one-dimensional linear potential (e.g., [T2"| 17]). 

Unlike in the quadratic problem, no special treatment is needed when both x and 
y equal 0. One simply has 

k 2 t 3 



12 



(107) 



Such a trajectory exists for any value of t (with a = kt = q(0), b = 0, according to (llOip 
and (11021) ). in marked contrast to the quadratic case. 

5.2. Type E 

Type E is simply type B attached to two type-A paths. Therefore, the action is 
S E {x, y, t) = S A ,i(x, y, t) + S B (x, y, t) + S A ,2{x, y, t) 



1/ 



+ 



x 



k 2 {t 2 -Uf 



4ti ' 4(t-t 2 ) 12 
In the potential region we have a parabolic trajectory of the form 

q{r) = —kr 2 + ar + c 



(108) 
(109) 



WKB Approximation to the Power Wall 



22 



that crosses the axis at r = t\ . We need a formula for t\ in terms of the path data, 
(x,y,t). Solving f ll09[) for the parameter a we obtain 

kt 2 - c 



ti 

Then the initial-velocity condition at t\ 



110) 



q(t 1 )=v 1 = -^ = -2kt 1 + a, (111) 



yields 



-y + 2kti 



2 



;ii4) 



c = y-ktt (112) 
t\ 

Substituting f )112p back into (jlU9j) and solving q(r) = 0, we find that the second axis- 
crossing time is 

kt? — y , 

sr- (113) 

On the other hand, since v 2 = —v\ , we have 
x _ y_ 
t-t 2 h' 

Combining (I113p and ( 11141) . we arrive at a quadratic equation for t\ , whose solutions 
are 

h{x, y, t) = ^—{t ± [t 2 + 4(x + y)/k} 1/2 }. (115) 
2(2/ + x) 

Existence requires 

«, 41a; + 7/1 

t 2 > 1 |L£l . 116 
k 

The qualitative reason for this condition is that short time lapses force q(ti) to be small 
for the type-B segment and large for the type-A segments, and these requirements come 
into conflict if t is too small. When the inequality ( II 16ft is strict, there are two allowed 
paths, one of low velocity that barely enters the potential region and one of high velocity 
that spends most of its time there. (For example, if x = y = —3kt 2 /32, one solution 
has t\ = 3i/8 and t 2 = 5t/8 and the other has t\ = t/8 and t 2 = 7t/8.) 

Substituting (11151) and f 1 1 1 3 j) into (I108p . we express the action just in terms of the 
(x, y, t) variables: 

S E (x, y, t) = ±(k 2 t 2 + Ak(x + y)f' 2 - - ^ . (117) 

It follows that 

A 2_ -yft _ 5 

2^ + 4(x + y)' A [kt 2 + A{x + y)] 2 ' 1 1 

Clearly the solution is not exact (A A ^ 0), and in fact the path has passed through a 
caustic, as shown by the negative sign in f 1 1 1 8 j) and the analysis in the next subsection. 
When (I116P is an equality, x is actually sitting on the caustic point. 
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Q 



Figure 8. Type E: (Linear Potential Case) The trajectory in the middle is not a 
sinusoidal path in this case, but rather a parabolic one. The maximum excursion 
occurs at time i m id , halfway between the times when the path crosses the axis. 



For use in the next subsection we note that if x = and also the caustic equality 
holds, then 



ti = t 2 = -^, t x =-t 2 , t? = - V -. (119) 
k 2 k K J 



5.3. Type C 



A path of type C starts at a point y in the free region and ends at a point x in the 
potential region. So the classical action is the sum of the action for the direct path and 
the action for a type-B path: 

Sc( X , »,«) = ■£ + - Xk(t ~ U) - k2{t ~ ' l)3 , (120) 

cv ' y ' ; 4ti 4(t-ti) 2 12 v ' 

where t is the final time and t% is the time when the particle passes through the origin. 

In the potential region the trajectory has the form (I109p and must satisfy 

x = q{t) = -kt 2 + at + c } - — = q(tx) = -2kt x + a. (121) 

t\ 

Hence 

a = 2kh - — , c = x + kt 2 - 2ktt 1 + — . (122) 
ti t\ 

Therefore, 

= q(t x ) = +kt 2 -y + x + kt 2 - 2ktt x + yt/t x , (123) 

which gives a cubic equation for t\ , 

kt\ - 2kt t\ + (kt 2 + x- y)t x + yt = 0. (124) 

With (j!24p the formula for c in (1 1 2 2 [) can be simplified, so that the final formula for the 
trajectory is 



q( T ) = -kr 2 + \ 2kt x - j- j r + (y - kt{). (125) 

We can rewrite equation (j!24p as 

t\ + a 2 t\ + a x t x + a = (126) 
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with 

a 2 = -2t, a 1 = ^—^- + t 2 ) ao = j- (127) 

A cubic equation always has three roots. Roots can be either all real, with multiplicity 
1, 2, or 3, or one real with the others complex conjugates (in the case of real coefficients). 
The nature of the roots is indicated by a polynomial discriminant [25] UHl [2U] [D P- 17] [23, 
pp. 103-105]. In terms of the auxiliary quantities 

3ai-a§ 9a 2 ai - 27a - 2a^ 
Q = g , R= 54 , (128) 

the cubic equation becomes 

r 3 + 3Qr - 2R = 0, r = t - § t, (129) 
and the discriminant is defined as 

D = i? 2 + Q 3 

= (3A;)" 3 [t 4 fc 2 x + t 2 k (2x 2 + hxy - \y 2 ) + (x - y) 3 ] . (130) 

If D > 0, one root is real and the other two are complex conjugates. If D = 0, all the 
roots are real with at least two equal. And if D < 0, then all roots are real and distinct. 
We note that when t = 0, 

(3&) 3 £> = (x - yf > (131) 

(since y > 0), and also that 

lim D(t) = +oo. (132) 

t— )-±oo 

So, as a function of t, either D{t) is everywhere positive or it dips negative for some 
interval of t before becoming positive again. In the latter case it must have a minimum. 
To analyze this behavior, we note that only positive t values are of interest, so we can 
examine the equation as a function of t 2 without loss of generality. Letting T = t 2 , we 

see 

(3A;) 3 ^ = 2k 2 xT + k (2x 2 + 5xy - \y 2 ) . (133) 
Setting (I133p to zero determines the positive critical point of D to be 



-8x 2 -20xy + y 2 , , 

t = V sr^- < 134 > 

Existence of this critical point requires 

C = -8a; 2 - 20sy + y 2 > ; (135) 

if (11351) is violated, D is always positive. When (11351) holds, substituting ( II 341) into the 
discriminant (11301) gives 

y(8x + y) 3 
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This expression can have any sign, since x > but y < 0. This shows that all possible 
cases of cubic roots can occur (unlike the related situation in [26] , where there are always 
3 real roots but only the middle one is physically relevant). 

Now we define s = 8x + y, as suggested by f)136p . and change variables from (y, x) 
to (s, x). Note that — oo < s < 8x. From (11351) we have C = s 2 — 36xs + 216a; 2 . Solving 
(11301) as a quadratic equation in T = t 2 , we obtain 

T = tk^ c± ^\ < 137 > 

where the new discriminant is 

P = 32x 3 + 12x 2 y 2 + fa;?/ 3 + ±y 4 

= -l xs3 + re 3 * = T6 s3 y- ( 138 ) 

If P is negative (0 < s < 8x), then there are no real roots for T and again D > for 
all T. If P > (s < 0), then there are two real roots, and they are positive because 
C > in this case; thus D < for some range of t. It can be seen that C(s) = has 
two roots, one slightly less than 8x and one in the irrelevant region s > 8x. Thus the 
results of this analysis can be listed: 

• s < =>- 3 real roots exist for some t. 

• < s < s_ w 7.6a; =>■ 1 real root for all t. 

• s_ < s < 8x =>- 1 real root for all t. 

• s>8x =3- y > (forbidden). 

The only distinction between the second and third cases is that in the second, D m[n is 
greater than 0, while in the third, the minimum does not exist at all. 

When there is only one real root, it is conveniently represented by the purely 
algebraic formula developed by Del Ferro, Tartaglia, Cardano and Bombelli (5j pp. 
310-317]: 

x = \t+ sJr+ ^R2 + Q3 + ?Jr- ^/WT^. (139) 

When there are three real roots, however, formula (I139p . although still correct, is almost 
useless in practice. A better formula is the trigonometric representation developed by 
Viete and Girard [5j p. 341]: 

ri (x, y,t) = ~t + 2 v /-g(x, j/ ,t)cos( 9( Y ,t) ) , (140) 

r 2 (x, y, t) = ~t + 2 V'Q(x,y,t)cos{ e ^ y f + 2n ) , (141) 

r 3 (x, y, t) = 2 -t + 2 v/-Q(x,y,t)cos( 9(x ^ 3 t)+47r ) , (142) 

where 
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Note that r 2 = r 3 when = 0, which is the same as the boundary D = between the 
two regimes. 

For some purposes, however, it is convenient to ignore the complicated explicit 
solutions (I139p and (I140p - (I143I) and work by implicit differentiation as in section 14.41 
We write 

/(^) = kt\ - 2kt t\ + (kt 2 + x - y)ti + yt (144) 

for the cubic polynomial in ( 1124ft . In calculating A 2 we shall encounter which 
involves the reciprocal of the derivative of the implicit equation /(ti) = with respect 
to t\ . A caustic may arise when this factor vanishes. As in the case of the quadratic 
potential, a caustic is signaled by a simultaneous solution of 

f(ti) = 0, f (t x ) = 0, (145) 

i.e., a place where the graph of / is tangent to the horizontal axis. But this is the same 
as saying that the cubic equation has a double root. So a caustic occurs at a point where 
the number of real roots is changing from one to three as a parameter of the problem 
varies. We define 

Y = f (t0 = 3H! 2 - Aktt x + kt 2 + x-y. (146) 
Differentiating (11241) with respect to x yields 



and thus 



Therefore, 



= t 1 +Y- 

ox 



dh ti 
dx Y 



, d 2 S dp y dti y 



(147) 



dxdy dx 2t 2 dx 2t x Y ' 1 } 



Then using (J84J), and using the cubic equation repeatedly to reduce the numerator to a 
quadratic in t\ , we obtain 

A A k 

r -24t 2 y - 9t(kt 2 + 2x - y)ti + (23kt 2 - 2lx + 2ly)t?]. (149) 



A Y A 

Finally, we note an interesting special case. In the equation D = 0, set x = (i.e., 
consider a caustic occurring right on the vertical axis). Then t 2 — t, and from (11301) we 
have 

= (3A;) 3 {-^f- V 3 
or 

y = - — ■ ( 15 °) 

But in this situation we also have 

«-«. = -£, (151) 
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because the final remark in section I5TT1 obviously generalizes to the statement that q{t\) 
is equal to k times the elapsed time of the excursion through the potential region. 
Substituting (I15ip into (11501) yields for y the equation 

y 2 + 2kt? + k\ 4 = 0, 

whose only solution is y = —kt\. By comparison with fl 1 5 [) it follows that t — 2t\ . 
Thereby we have reproduced the relations (11191) obtained by approaching the axis caustic 
from the other side. 



5.4- Type D 



The part of a type D trajectory in the potential region is determined by the initial 
conditions g(0) = y and q{ti) = 0. Thus 

ktl-y 



-kr + ar + y, a 



ti 



(152) 



Also, the final veolcity must be the same as the velocity when the particle crosses the 
vertical time axis: 

x 



(t - h) 



q(ti 



-2kti + a. 



Hence 



x 



+ 2/cti . 



(t-ti) 

The Lagrangian of the particle in the linear potential region is 

L(x, y , t ) = *M-^ + ^_ 2k\r + 2 -^L + 2k\\ 



so the action for the particle in that region is 



S B (x,y,t) 



L(q, q, r) dT = -—kt x — + 



2 

4tT 



The total action is the sum of the actions for the two path segments: 



S D (x, y, t) = S A (x, y, t) + S B (x, y, t) 



4(t-ti) 12 1 2 



2 

4*7 



The final momentum of the classical particle is given by p = — |f , and hence 



i x 

V = —V = —, r . 

F 2 2(i-ii) 



Taking the partial derivative of momentum with respect to x gives 



or 



A 2 



A 2 



d 2 S dp d 



x 



dydx dy dy 12 (t — ti) 
9ti 



x 



2{t-t 1 fdy 



(153) 



(154) 



(155) 



(156) 



(157) 



(158) 



(159) 
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The partial derivative of t\ with respect to x can be found from 



q{h) = 
which simplifies to 



-kt\ + 



x 



(t - h) 



+ 2Mi 



h + y, 



= kt\{t - t x ) + xh + y(t - h). 

Implicit differentiation then yields 

dti ti dti t — ti 



dx X 



dy 



X 



where 



Thus 



X = -x + y - {2kt)t x + 3kt 2 . 



A 2 



x 



2(t-t 1 )X ' 

Finally, after a long Mathematica calculation we obtain 
AA 1 



(160) 

(161) 

(162) 

(163) 
(164) 



A x 2 (t-t0 2 X 4 

[t 2 y(-Ak 2 t 4 y - A7x 2 y + 76xy 2 - 20y 3 + 6kt 2 x 2 - 20kt 2 xy + A0kt 2 y 2 ) 
+ 2ty(-21x 2 + 12k 2 t A y + 69x 2 y - 60xy 2 + 12y 3 - 22kt 2 x 2 

+ 60kt 2 xy - A0kt 2 y 2 % 
+ (5x 4 + 10x 3 y + 38kt 2 x 2 - 39x 2 y 2 - 84kt 2 xy 2 + 28xy 3 

- 20k 2 t 4 y 2 + A0kt 2 y 3 - 4y*)t?]. (165) 



6. Conclusion and Outlook 

The WKB, or semiclassical, approximation to the propagator of a quantum system in 
more than one dimension has a well-developed literature. Here we have extended it in 
two directions. 

In section [2] we investigated the convergence of the semiclassical series and proved, 
under rather strong assumptions, that it converges, thereby giving a construction of 
the exact quantum propagator as a sum or integral over paths (classical except for 
scatterings off the potential at finitely many points). 

In later sections we constructed the leading term in the series for two instances 
of a special class of potentials, the soft walls of integral power growth. Such detailed 
applications to particular systems are rather rare in the literature, a fact that becomes 
less surprising once one discovers how complicated the calculations can be even for 
seemingly simple examples. The first step in such a venture is to find all the classical 
paths between given starting and ending points in space-time; the paths may not be 
unique (and may not exist), and they fall into several qualitatively different classes. For 
the particular cases of a = 2 (harmonic potential inside the wall) and a = 1 (linear 
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potential inside the wall) we studied all the paths in detail, constructing their associated 
semiclassical apparatus of action, amplitude, and residual function (proportional to the 
Laplacian of the amplitude). There are noteworthy differences between the harmonic 
and linear cases; generally speaking, the former is simpler because of the fixed period 
of harmonic oscillations. 

For a solution of the time- dependent Schrodinger equation, a caustic is not the same 
thing as a breakdown of the approximation; as the propagator for the harmonic oscillator 
shows, a singular function may be an exact solution. For the power walls, however, it 
turns out that "bad" caustics generically occur; that is, there is a singularity in the 
residual function as well as the amplitude. Thus the theorem from section [2] ceases to 
apply after a short initial time period. 

We have left many things undone that are suitable for future work. First, there 
is the possibility of "getting over the caustic" by constructing a first-order WKB 
approximation in momentum space. We hope that using this improved approximation 
in the higher-order WKB construction will yield a globally convergent series. 

Of course, we are far from calculating any such series even in the caustic-free region. 
Our solutions for paths of Types C and D are ready to be concatenated to create the 
first-order term in the series by integration over a scattering point inside the wall. 
Evaluating such integrals in practice may prove challenging. 

Another problem that has been skipped over here is the diffraction from the mild 
singularity in our potential at the origin. Wherever a potential fails to be of class 
C°°, the semiclassical expansion fails to be valid beyond a certain order. A correct 
approximation in the WKB spirit must then take into account exactly the scattering 
from the point in question. The resulting correction is small compared to the leading 
semiclassical approximation but will be more significant than some higher-order terms 
in the expansion fl28|) . For the power wall with integer exponent a, the discontinuity 
appears in the derivative of order a. For a = 1 it is quite significant, and its effects show 
up in the calculation of quantum vacuum energy in that potential background [TH H] . 

The physical problem motivating [I] and [19] is vacuum energy in relativistic 
quantum field theory, where the natural integral kernel to study is not the Schrodinger 
propagator but the Wightman or Feynman wave kernel (or its continuation to imaginary 
time, the cylinder or Poisson kernel). The latter is not so amenable to WKB techniques, 
and an asymptotic expansion for one kernel does not automatically yield one for the 
other. However, both kernels encode spectral information about the same elliptic 
operator acting in the spatial coordinates, and we continue to investigate how to use 
the nonrelativistic semiclassical construction in the study of the relativistic theory. 
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